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ABSTRACT 

Using a simple two-color selection based on g-, z-, and X-band photometry, we pick out 1609 
star-forming galaxies (sgzKs) and 422 passively evolving galaxies (pgzKs) at z ~ 2 from a i^-band- 
selected sample (Kab < 22.0) in an area of ~ 0.44 deg 2 of the All-wavelength Extended Groth Strip 
International Survey. The number counts of pgzKs in our sample turn over at Kab ~ 21.0, and both 
the number of faint and bright objects (including sgzKs and pgzKs) exceed the predictions of a recent 
semi-analytic model of galaxy formation, a more successful model is need to explain this diversity. 
We also find that the star formation rate (SFR) and specific SFR (sSFR) of sgzKs increases with 
rcdshift at all masses, implying that star-forming galaxies were much more active on average in the 
past. Moreover, the sSFR of massive galaxies is lower at all redshifts, suggesting that star formation 
contributes more to the mass growth of low-mass galaxies than to high-mass galaxies. From Hubble 
Space Telescope Wide Field Camera 3 near-infrared imaging data, we find that morphologies of z ~ 2 
galaxies not only have diffuse structures with lower G and higher Mao values, but also have single- 
object morphologies (higher G and lower M20), implying that there are morphological variety and 
different formation process for these galaxies at z ~ 2. Finally, we also study the fraction of active 
galactic nuclei (AGNs) in the gzKs, 82 of 828 gzKs with four IRAC bands can be classified as AGNs 
(~ 10%). Most of these AGN candidates have Ln.5-10 kcV > 10 41 ergs -1 . 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: formation - galaxies: 
high-redshift - galaxies: photometry 



1. INTRODUCTION 

Understanding when and how the most massive galax- 
ies in the universe formed is one of the most outstand- 
ing problems in cosmology and galaxy formation (Con- 
selice et al. 2007). A number of observations suggest 
that the era of z ~ 2 is important in galaxy forma- 
tion and evolution for various reasons: the cosmic star 
formation rate (SFR) density begins to drop at z ~ 2 
from a flat plateau at higher redshifts; the morpholog- 
ical type mix of field galaxies changes remarkably at 
z ~ 2; the number density of QSOs has a peak at z ~ 2; 
~ 50% — 70% of the stellar mass assembly of galaxies 
took place in the redshift range 1 < z < 3 (Dickinson 
et al. 2003; Fontana et al. 2003; Steidel et al. 2004; Kong 
et al. 2006; Richards et al. 2006; Arnouts et al. 2007; 
Pozzetti et al. 2007; Noeske et al. 2007). However, in the 
redshift range 1.4 < z < 3.0, identifiable spectral fea- 
tures such as Ca II H&K and the 4000 A break (for pas- 
sive galaxies) and [O II]A3727, [O III]A5007, Ha, and H/3 
(for star-forming ones) move out of the optical bands, 
thus near-infrared imaging and spectroscopy become es- 
sential. On the other hand, near-infrared galaxy sam- 
ples offer several advantages compared to purely optical 
selections (Cowie et al. 1994). They allow us to select 
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galaxies at z > 1 in the rest-frame optical, correspond 
more closely to a stellar-mass-selected sample, and are 
less prone to dust extinction. 

There are many ways of using optical and infrared 
(IR) colors to select the sample of galaxies at z ~ 2. A 
sample of IR- luminous dust-obscured galaxies (DOGs), 
with a median iyL L ,(8fim) ~ 4 x 10 11 Lq, selected to 
have very red R — [24] color, was recently presented 
by Dey et al. (2008). From spectroscopic observations 
DOGs were found to have a tight redshift distribution 
around z ~ 2, very similar to that of the submillimeter- 
selected galaxies (Chapman et al. 2005). Another differ- 
ent sample, based on near-infrared color selection, is that 
distant red galaxies with ( J— K ) VO ga > 2.3 (DRGs; Franx 
et al. 2003; van Dokkum et al. 2004). These sources are 
redder and strong clustering and are believed to be 
more massive than 10 11 M Q (Labbe et al. 2005; Papovich 
et al.2006; van Dokkum et al.2006; Quadri et al.2008). 
A third sample uses Spitzer/IRAC color to select galaxies 
at z > 1.4. In the redshift range 1.4 < z < 2.7, the four 
IRAC bands probe the rest-frame near-infrared bands 
where galaxy spectral energy distributions (SEDs) have 
similar shape, thus the IRAC colors are very robust in 
determining rcdshift in this range (Farrah et al. 2008; De- 
sai et al. 2009; Huang et al. 2009; Fadda et al. 2010; Fi- 
olet et al. 2010; Fang et al. 2011). Most sources of these 
samples are starburst-dominated ultraluminous infrared 
galaxies (ULIRGs; Lir > 10 12 Lq), with an apparent 1.6 
/jm stellar bump in the IRAC channels. 

An independent way to select galaxies at z ~ 2 is to use 
BzK colors. Daddi et al. (2004) introduced a new photo- 
metric technique, BzK = (z — K)ab — (B — z)ab = —0.2 
and (z — K)ab = 2.5, for obtaining a virtually com- 
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plete sample of star- forming (BzK > —0.2) and passively 
evolving galaxies (BzK < —0.2 and (z — K)ab > 2.5) at 
1.4 < z < 2.5. These criteria arc reddening independent 
for star-forming galaxies in the selected redshift range, 
thus can be used to select the reddest dust-extinguished 
galaxies. This should allow a relatively unbiased se- 
lection of star-forming (sBzKs) and passively evolving 
galaxies (pBzKs) at z ~ 2 within the magnitude limit 
of the sample studied. Based on BRIzJK photome- 
try from Subaru and New Technology Telescope (NTT) 
over two separate fields, using the BzK technique, Kong 
ct al. (2006) have obtained a sample of ~ 500 sBzKs and 
~ 160 pBzKs, which were identified over an area of ~ 920 
arcmin 2 to Kab — 20.8, of which 320 arcmin 2 are com- 
plete to Kab = 21.8. They found that the log of the 
number counts of pBzKs flattens out by Kab ~ 20.8. 
Similar conclusions were also found by Lane et al. (2007) 
and McCracken et al. (2010), with much larger samples. 

Investigations of the physical properties (i.e., stellar 
mass, luminosity, size, clustering, SFR, morphology, etc.) 
of various types of galaxies from high to low redshift 
are required for better understanding the history of the 
galaxy formation and evolution. Dunne et al. (2009) 
studied the star formation history of ^-selected galax- 
ies in the UKIDSS Ultra-Deep Survey. They found that 
the specific star formation rate (sSFR) for if-selected 
sources rises strongly with redshift at all stellar masses. 
In the same field, Williams et al. (2009, 2010) investi- 
gated the physical properties (including clustering, size, 
stellar mass, surface density, SFR, color, etc.) of qui- 
escent and star-forming galaxy populations to z =0 
2 with pure photometric data, employing a novel rest- 
frame U — V versus V — J technique (Wuyts et al. 2007; 
Wolf et al. 2009; Balogh et al. 2009; Whitaker et al. 2010; 
Patel et al. 2011). As shown in Figure 8 of Williams 
ct al. (2009) and Figure 2 of Brammer et al. (2011), the 
red quiescent and star-forming galaxies are found to oc- 
cupy two distinct population in the rest- frame U — V ver- 
sus V — J color space. This bimodal behavior is still seen 
up to z ~ 2. At all redshifts, furthermore, massive quies- 
cent galaxies occupy the extreme high end of the surface 
density distribution and a tight mass-size correlation, 
while star-forming show a broad range of both densi- 
ties and sizes. On the other hand, the sizes and surface 
densities of massive quiescent and star-forming galaxies 
evolve as simple power-law behavior in (1 + z), with 
more massive galaxies exhibiting faster evolution. In 
the meantime, Brammer et al. (2009, 2011), van Dokkum 
ct al. (2010), and Whitaker et al. (2011) provided fur- 
ther evidence for the results above, using the sample 
from the NEWFIRM Medium Band Survey (NMBS; van 
Dokkum et al. 2009). More studies for the 7-sT-selected 
galaxies have been performed up to z = 5 (Cimatti 
et al. 2002; Abraham et al. 2004; Drory et al. 2005, 2009; 
Fontana et al. 2006; Pozzetti et al. 2007; Perez-Gonzalez 
et al.2008; Furusawa et al.2011). 

Morphologies are also essential in studying galaxy mass 
assembly history and evolution. For the morphologies of 
galaxies at z ~ 2, the traditional Hubble sequence of 
regular spirals and elliptical galaxies has not settled into 
place by z ~ 2, and a much higher frequency of diffuse, 
clumpy, and irregular structures is observed among star- 
forming systems (Lotz et al. 2006; Szomoru et al.2011). 
With the advent of Hubble Space Telescope HST/Wide 



Field Camera 3 (WFC3), it is now possible to obtain 
a complete rest-frame optical morphological census for 
massive galaxies at z ~ 2. 

Understanding the nature of active galactic nuclei 
(AGNs) and the galaxies that host them is important 
for such diverse goals as pinpointing the sources of the 
cosmic X-ray and IR backgrounds and deriving the star 
formation history of the universe. Multi-wavelength sur- 
veys are particularly important for the study of AGNs 
because their appearance in different wavelength regimes 
can be quite different. In this work, we will discuss the 
fraction of AGNs in our sample by using the IRAC band 
color criteria of Stern et al. (2005) and mid-IR spectral 
index (Barmby et al. 2006; Park et al. 2010). 

To study the physical properties of galaxies in the 
redshift range 1.4 £ z £ 2.5, in this paper, we select a 
sample of passively evolving and star-forming galaxies 
to Kab < 22 in the All-wavelength Extended Groth 
Strip International Survey (AEGIS) data set. The pa- 
per is organized as follows. We describe the multi-band 
spatial- and ground-based observations of the AEGIS 
field, and introduce data reduction, photometric red- 
shifts, and method for galaxy sample selection in Sec- 
tion 2. Section 3 presents source number counts (sgzKs 
and pgzKs). Section 4 presents the physical properties 
of sgzKs (including the estimate SFR and stellar mass 
of star-forming galaxies and the stellar mass versus SFR 
relation). Section 5 describes the morphologies of both 
sgzKs and pgzKs in the AEGIS. Section 6 discusses the 
fraction of AGNs in our galaxy sample. Finally, a brief 
summary is presented in Section 7. All magnitudes and 
colors are in the AB system unless stated otherwise. 
Throughout the paper we adopt the following cosmol- 
ogy: h = H [km s- 1 Mpc-^/lOO = 0.71, Q A = 0.73, 
VL M = 0.27. 

2. DATA AND SAMPLE SELECTION 

The AEGIS is a collaborative effort to obtain both 
deep imaging covering all major wavebands from X-ray 
to radio and optical spectroscopy over a large area of sky 
(0.5 — 1.0 deg 2 ) with the aim of studying the panchro- 
matic properties of galaxies over the last half of the 
Hubble time. AEGIS is targeted on a special area of 
the sky, called the Extended Groth Strip (EGS; centered 
at a(J2000)= 14 h 19 m 12 s , (5(J2000)= 52°43'42"), which 
has low-extinction, low-Galactic, and zodiacal infrared 
emission, has good schedulability by space-based obser- 
vations, and has therefore attracted a wide range of deep 
observations at essentially every accessible wavelength 
over this comparatively wide field (Davis et al. 2007). 

2.1. Data Observations 

The EGS is a public deep multi- wavelength survey; this 
field provides an unique combination of area and depth 
at almost every waveband observable. In this paper, 
Chandra/ ACIS X-ray, CFHT/MegaCam Legacy Survey 
optical, CFHT/CFH12K optical, Palomar/WIRC near- 
infrared, Spitzer/IRAC mid-infrared, and Spitzer/MlPS 
far-infrared imaging data or catalog of the EGS is used. 
Davis et al. (2007) provided more details on the data re- 
duction and its performance. 

For our analysis, specially, the Chandra/ ACIS DR2 
catalog presented by Laird et al. (2009) is used. The X- 
ray data reduction was performed using the CIAO data 
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analysis software version 3.3. Basic data reduction pro- 
ceeded in a manner similar to that described by Nan- 
dra et al. (2005). The vast majority of X-ray sources de- 
tected in the EGS field are AGNs. In the meantime, 
the EGS is also one of four 1.0 deg 2 fields covered by the 
CFHTLS deep survey (labeled CFHTLS-D3). The data 
covers the observed wavelength range 3500A< A <9400A 
in the u*, g', r', i 1 , z' filters (hereafter labeled u, g, 
r, i, z). The limiting AB magnitudes is u ~ 27.0, 
g ~ 28.3, r ~ 27.5, i ~ 27.0, z ~ 26.4 in the current data 
set (corresponding to the magnitude limit at 5er point- 
source detection). Near-infrared data catalog in J and 
K bands are available over ~ 0.67 deg 2 (if-band) with 
the 5cr magnitudes limits of Jab ~ 23.9 and A"ab ~ 22.5 
(Bundy et al. 2006). Deep observations of the EGS with 
Spzfaer/Infrared Array Camera (IRAC) cover an area 
of ~ 0.33 deg 2 to a 50% completeness limit of 1.5 /xJy 
at 3.6 /xm (Barmby et al. 2008). The catalog comprises 
57,434 objects detected at 3.6 ^m, with 84%, 28%, and 
24% also detected at 4.5 /mi, 5.8 /j,m, 8.0 lira. As for far- 
infrared data sets, we also use the 5pzteer/Multiband 
Imaging Photometer for Spitzer (MIPS) 24 /im imaging 
data with the 5cr flux limit of ~ 77/xJy over 0.534 deg 2 . 
Sources are identified and photometry extracted with 
point-spread function (PSF) fitting using the DA0PH0T 
software (Stetson 1987). 

A catalog of objects with redshifts known from spec- 
troscopic observations is very important to the construc- 
tion of a photometric redshift code. For this reason, 
we use the DEEP2 (Deep Extragalactic Evolutionary 
Probe) spectroscopic data obtained with the Deep Imag- 
ing Multi-Object Spectrograph boarded on the Keck 10 
m telescope. The DEEP2-DR3 redshift catalog struc- 
tures contain 24 tags for each object with 50,386 en- 
tries in the entire catalog. For the EGS field (DEEP2- 
ficldl), we only consider objects with redshift quality 
greater than or equal to 3. To specifically study the 
accuracy of the photometric redshift, we also adopt spec- 
troscopic redshifts drawn from the Lyman break galaxy 
(LBG; z ~ 3) catalog of Steidel et al. (2003) and from 
a small sample of ULIRGs (z ~ 2) presented in Huang 
et al. (2009) and Fang et al. (2011). 

The NMBS employs a new technique of using 
medium-bandwidth near-infrared filters to sample the 
Balmer/4000A break from 1.5 < z < 3.5 at a higher 
resolution than the standard broadband near-infrared fil- 
ters (van Dokkum et al. 2009), thereby improving the ac- 
curacy of photometric redshifts. A custom set of five 
medium-bandwidth filters in the wavelength range of 
1 — 1.8 /tm were fabricated for the NEWFIRM camera 
on the Mayall 4 m telescope on Kitt Peak for the NMBS 
(see Figure 1 of Whitaker et al. 2011). The survey tar- 
gets two fields within the COSMOS (Scoville et al. 2007) 
and AEGIS (Davis et al. 2007) surveys, chosen to take 
advantage of the wealth of publicly available ancillary 
data over a broad wavelength range. The more details 
of the observations, data reduction, and photometry for 
the survey can be found in Whitaker et al. (2011). The 
NMBS catalogs contain ~ 13,000 galaxies at z > 1.5 
with accurate photometric redshifts and rest-frame col- 
ors. In this paper, we will use their photometric redshifts 
in the AEGIS field. 



2.2. Data Reduction 

Source extraction for the CFHTLS optical images is 
performed with SExtractor (Bertin & Arnouts 1996) in 
the dual image mode (i band is used as the reference 
image); a 2"diameter aperture is used for aperture mag- 
nitudes. The edges of the resulting final images, where 
the signal to noise of CFHTLS /MegaCam images is very 
low, are trimmed. We finally cross-correlate the i-band 
catalog (that included u, <?, r, and z photometry) with 
the -ftf-band catalog, with a matching radius 2". The fi- 
nal overlapping area of all our data sets is ~ 0.44 deg 2 , 
and 16,228 objects arc included in our final catalog 
(K AB < 22 mag). Comparing to the NMBS in the 
AEGIS field (~ 0.21 deg 2 ), we have larger area and more 
faint galaxies (see Section 2.4). Stellar object are isolated 
with the color criterion (z-A^)ab < 0.4(5 — z)ab~ 0.45, 
similar as Daddi et al. (2004). A total of 1847 objects 
are classified as stars, and the number of the final galaxy 
sample is 14,381. All magnitudes are corrected for Galac- 
tic extinction. 

2.3. Photometric Redshift 

An SED fitting technique based on an updated version 
of the HyperZ code (Bolzonella et al. 2000) is used to 
calculate the photometric redshifts (hereafter z p ) for all 
the Kab < 22 galaxies in the EGS. A set of templates 
generated by Brammer & van Dokkum (2007) for their 
EAZY photometric redshift code are used. For each ob- 
ject, the SED derived from the observed magnitudes is 
compared to each template spectrum in turn, and the 
weighted mean redshift, computed in the confidence in- 
tervals at 99% around the main solution, from HyperZ, 
is calculated. 

To check the accuracy of the photometric redshift 
we derived, we first compare the z p with the spectro- 
scopic redshifts. For the 14,381 galaxies in our sam- 
ple, 2940 of them have spectroscopic redshifts (z s ) from 
the DEEP2-DR3 redshift catalog, which have a confi- 
dence level greater or equal to 97%. In the meantime, 
spectroscopic redshifts are also available for 35 LBGs at 
z ~ 3 and 14 ULIRGs at z ~ 2 within the EGS field 
(14 ULIRGs presented in Huang et al. 2009 and Fang 
et al. 2011, 35 LBGs drawn from the catalog of Steidel 
et al. 2003). In Figure Ufa), we show a comparison of 
the photometric redshift with the spectroscopic redshift 
for those 2989 galaxies. Filled triangles and stars corre- 
spond to ULIRGs and LBGs with spectroscopic redshifts, 
respectively. From this figure, we find that our photomet- 
ric redshifts are in good agreement with the spectroscopic 
redshifts, with an average (z p — z s )/(l + z s ) = —0.014. 
Only 128 of 2989 galaxies have (z p -z s )/(l+z s ) > 0.1; the 
percentage of catastrophic failures is about 4.3%. The 
normalized median absolute deviation, ctnmad, f° r °ur 
galaxy sample is onmad = 0.032. 

On the other hand, Whitaker et al. (2011) described 
the full details of the observations, data reduction, and 
photometry from the NMBS in the AEGIS field (~ 
0.21 deg 2 ), also presented a public if-selected photomet- 
ric catalog along with accurate photometric redshifts and 
rest-frame colors. For this catalog, we will use their pho- 
tometric redshifts (znmbs)- We finally cross-correlate 
our catalog with the NMBS catalog in the AEGIS field, 
with a matching radius 1". A total of 5533 galaxies 
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Z s Z NMBS 

Fig. 1. — (a) Comparison between spectroscopic redshift (z s ) and photometric redshift (z p ). The dot-dashed lines are for z p = z s ± 
0.10(l + z a ). The number of galaxies with spectroscopic redshift and the fraction of catastrophic failures ((z p — 2 S )/(1 + z s ) > 0.1) are listed 
in the top-left corner. The inset shows the Az = z p — z s distribution, and the la dispersion is 0.032. The dashed line is z p = z s . Filled 
triangles and stars represent ULIRGs at z ~ 2 and LBGs at z ~ 3 respectively. Filled circles are drawn from the DEEP2-DR3 redshift 
catalog, (b) Comparison between znmbs an d zour- The znmbs are drawn from a public K-selected photometric catalog with photometric 
redshifts in the NMBS/AEGIS field (Whitaker et al. 2011). A total of 5533 galaxies are located in the overlapping area between NMBS 
and AEGIS. 



are located in the overlapping area between NMBS and 
AEGIS. In FigureHfb), our photometric redshifts (zqur) 
are in good agreement with NMBS photometric redshifts 
( znmbs) in the AEGIS field, although there are some 
large dispersion. 

2.4. z ~ 2 Galaxy Selection 

To select a catalog of objects at z ~ 2, we use the 
BzK selection technique which was introduced by Daddi 
et al. (2004). This technique has been adopted and tested 
in several subsequent studies (e.g., Kong et al. 2006; 
Lane et al. 2007; Blanc et al. 2008; Dunne et al. 2009; 
McCracken et al. 2010; Onodera et al. 2010; Yoshikawa 
et al. 2010). To make the comparison possible with pre- 
vious studies, we wanted our photometric selection crite- 
rion to match as closely as possible as the original BzK 
selection introduced in Daddi et al. (2004). Figure Ufa) 
shows the transmission profiles of the VLT-BzK filters, 
which were used in Daddi et al. (2004), the transmission 
profiles of the CFHT-gz, and the Palomar-if filters of 
AEGIS, which are used by this paper. From this figure, 
we find that the shapes of those filters are different. To 
account for the differences, we carefully compare the stel- 
lar sequence in the AEGIS to that of Daddi et al. (2004), 
with the similar procedure outlined in Kong et al. (2006), 
and apply small color terms to B — z and g — z, z — K 
in AEGIS and in Daddi et al. (2004), in order to obtain 
a fully consistent match. 

From this analysis, we obtain our new color criteria 
g — z and z — K to select z ~ 2 galaxies: gzK — (z — 
K)ab — l-4(<? — z)ab > 0.2 can be used to select z > 1.4 
star-forming galaxies (hereafter named sgzKs), gzK < 
0.2 and (z — K)ab > 2.7 can be used to isolate z > 1.4 
passively evolving galaxies (hereafter named pgzKs) . To 



further elucidate the validity of these gzK criteria, we 
generate a set of templates, with different star formation 
history, by Bruzual & Chariot (2003) stellar population 
synthesis models, and find that galaxies at 1.4 ;$ z ^ 2.5 
with ongoing star formation indeed lie in the gzK > 0.2 
region as expected, passively evolving galaxies lie in the 
gzK < 0.2 and (z — K )ab > 2.7 region of the g — z and 
z — K diagram. 

FigureEj^b) shows the gzK color diagram of -ftT-selected 
objects in the AEGIS. Using the color criterion gzK — 
(z - K)as - l-4(.g - z)ab > 0.2 (the black solid line), 
1609 galaxies with Kab < 22 are selected as sgzKs (blue 
squares), which occupy a narrow range on the left of 
the solid line in Figure EJb). Using gzK < 0.2 and 
(z — K)ab > 2.7 (dot-dashed line), 422 objects are se- 
lected as candidate pgzKs (red circles), which lie in the 
top-right part of Figure^b). The sample of gzKs (pgzKs 
and sgzKs) with open circles from the overlapping area 
between NMBS and AEGIS. These gzKs mostly lie along 
the boundary with non-gzK objects, comparing to the 
photometric redshift distribution for gzKs in our sam- 
ple, their photometric redshifts have lower values. The 
surface density of Kab < 22 sgzKs + pgzKs is ~ 1.28 
arcmin -2 , which is in good agreement with those pre- 
sented in Daddi et al. (2004), about ~ 1.1 arcmin -2 , and 
Kong et al. (2006), about ~ 1.4 arcmin -2 . 

In Figure [3ja) , we show the histogram of the photo- 
metric redshift for all iiT-selected galaxies in our sample. 
To check how efficient is our gzK color criteria in singling 
out z > 1.4 galaxies, the cross-hatched region shows all 
non-gzK galaxies in our sample. From this figure, we 
find that only 4% of galaxies with photometric redshift 
z > 1.4 cannot be selected by the gzK color criteria, 
which means the gzK criteria are quite efficient in sin- 
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Fig. 2. — (a) Filter transmission profiles for the AEGIS and those used by Daddi et al. (2004). These profiles are normalized to a 
maximum throughput of 1. The solid lines correspond to the B-, z-, and _R" s -band filter curves used in Daddi et al. (2004), respectively; the 
dashed lines correspond to the </-, z-, and i^-band curves used in the AEGIS, (b) Two-color (z — K)a_b vs. (g — z)ab diagram for the galaxies 
in the EGS field. Galaxies at high redshifts are highlighted. The diagonal solid line defines the region gzK = (z — K)ab — 1.4(g — z)ab > 0.2 
that is efficient to isolate z > 1.4 star-forming galaxies (sgzKs). The horizontal dot-dashed line further defines the region (z — i^)AB > 2.7 
that contains old galaxies at z > 1.4 (pgzKs). The dashed line defines the region occupied by stellar objects. Filled stars show stellar 
objects in our fields with (z — ~ 0.4(g — z)ab < —0.45; filled squares represent sgzKs; filled circles represent pgzKs. The sample of 

gzKs (pgzKs and sgzKs) with open circles from the overlapping area between NMBS and AEGIS. These gzKs mostly lie along the boundary 
with non-gzK objects; comparing to the photometric redshift distribution for gzKs in our sample, their photometric redshifts have lower 
values. 

number counts of pBzKs, in general, are in good agree- 
ment with those presented in the literature for z > 1.4 
passively evolving galaxies. At bright magnitudes (e.g., 
Kab < 21.0) our counts are in very good agreement 
with the pBzKs counts from McCracken et al. (2010). 
However, the number counts of ours and those of Mc- 
Cracken et al. (2010) are slightly below those from Blanc 
et al. (2008) at faint magnitudes. Combining literature 
data with the present work, we can examine the global 
shape of the z > 1.4 passively evolving galaxy counts 
over a very wide dynamic range. This strongly suggests 
a break feature in the slope at Kab < 21.0, which was 
found in Kong et al. (2006) already. 

The dot-dashed lines in Figure HJa) and (b) show the 
number counts from the semi-analytic model of Kitzbich- 
ler & White (2007). The number counts of quiescent 
galaxies and star-forming galaxies were derived by Mc- 
Cracken et al. (2010), following the approach of Daddi 
et al. (2007). From Figure HI we find the models over- 
predict the number of faint galaxies, but underpredict 
the number of bright galaxies in all two plots. Given 
the narrow redshift range of our z > 1.4 galaxy popu- 
lation, apparent K magnitude is a good proxy for ab- 
solute K magnitude, which can itself be directly related 
to underlying stellar mass (Daddi et al. 2004) . This im- 
plies that these models predict too many small low-mass 
star-forming and passively evolving galaxies, and too 
few large high-mass star-forming and passively evolving 
galaxies at z > 1.4. Similar conclusions were drawn 
by Fontana et al. (2006) and McCracken et al. (2010). 
Therefore, models incorporating AGN feedback similar 
to Kitzbichler & White undcrpredicted the number of 



gling out high-redshift galaxies. Figure E^b) shows the 
distribution of the photometric redshifts of sgzKs, with 
a median value z ~ 1.8. Most of them have redshifts at 
the range 1.4 < z < 2.5, only 8.5% of them with z < 1.4 
and 8.3% with z > 2.5. We show the redshift distribu- 
tion of pgzKs in Figure [3Jc) ; it clearly shows that there 
is a considerable fraction (12.8%) of objects at low red- 
shift (z < 1.0). The median redshift value of pgzKs is 
~ 1.5, which is similar to that (z ~ 1.4) of McCracken 
et al. (2010) . 

3. NUMBER COUNTS OF GZKS 

Figure |4] shows the differential number counts of sgzKs 
(left panel) and pgzKs(right panel). In both pan- 
els of Figure. 21 results from the Deep3a-F (asterisks) 
and Daddi-F (open diamonds) from Kong et al. (2006), 
UKIDSS (open squares, Lane et al. 2007; crosses, Hartley 
et al.2008), MUSYC (open triangles, Blanc et al.2008), 
and COSMOS (open circles, McCracken et al. 2010) are 
also plotted for comparison. Figure IHa) shows the num- 
ber counts of sgzKs, no correction for incompleteness is 
applied to our data. As shown in the figure, our num- 
ber counts of sgzKs in the EGS field are in good agree- 
ment with the number counts of z ~ 2 star-forming 
galaxies in Kong et al. (2006), Blanc et al. (2008), and 
McCracken et al. (2010). However, the counts presented 
by Lane et al. (2007) are offset compared to our counts, 
which is also found by McCracken et al. (2010), and the 
cause for this discrepancy may arises from an incorrect 
transformation to Daddi et al. 's (2004) filter set in Lane 
et al. (2007). 

Figure HJb) shows the number counts of pgzKs. Our 
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Fig. 3. — Photometric redshift distributions for the if -selected galaxies in our sample. The left panel shows the redshift distributions 
for all galaxies and all non-gzKs galaxies (cross-hatched). The top-right panel shows the redshift distributions of our sgzKs sample (1609 
galaxies), the bottom-right panel shows the redshift distributions of pgzKs (422 galaxies). 



high-mass galaxies, however predicted too many low- 
mass galaxies. 

4. PHYSICAL PROPERTIES OF sgzKs 

In this section, physical properties of sgzKs, such as 
stellar mass (M»), SFR, and sSFR (SFR/M*), are de- 
rived on the basis of the present photometric data. 

4.1. Stellar Mass Estimates 

In order to estimate the stellar mass (M») of sgzKs, 
we perform SED fitting of the multi-band photometry 
as described above with population synthesis models. 
The SED templates are generated with the SPS pack- 
age developed by Bruzual & Chariot (2003). We assume 
an universal initial mass function (IMF) from Chabrier 
(2003). Masses derived assuming this IMF can be con- 
verted to Salpeter by adding 0.3 dex (Bundy et al. 2005). 
A grid of model templates spanning a range of star for- 
mation histories (parameterized as an exponential), ages, 
metallicities and dust content are used. The stellar mass- 
to-luminosity (M/L) ratios from the best-fit templates 
and the photometric redshifts (spectroscopic redshifts are 
adopted if available) from Section 2.3 are used to calcu- 
late the stellar mass. The histogram for the stellar mass 
of the sgzKs is shown in Figure [SJa). About 50% of the 
sgzKs in our sample have M* > 10 n Af Q , and the median 
stellar mass of sgzKs is ~ 8.8 x 10 10 M Q . 

4.2. Star Formation Rates 

The ultraviolet continuum can be used to estimate 
the luminosity from young stars and hence the SFR. 
We derive the SFR of each sgzKs based on its rest- 
frame 2800A luminosity, which is calculated from the 
same best-fitting SED template as used in the estima- 
tion of the photometric redshift. Since our multi-band 
photometric data completely cover the rest-frame 2800A 



for the redshift range (1.0 ;$ z ~ 3.0) of sgzKs, the es- 
timated rest-frame 2800A luminosity is a robust quan- 
tity, especially in the case where a galaxy is detected in 
both bands which straddle the rest-frame 2800A. The 
color excess, E(B — V), from the SED fitting is used 
to derive the dust-corrected UV luminosity. Then, we 
convert the dust-corrected rest-frame 2800A luminosi- 
ties into SFRs using the calibration by Kennicutt (1998): 
SFR (M Q yr- 1 ) = 1.4 x lO" 28 L v (erg s" 1 Hz" 1 ). 

The SFR histogram of the sgzKs is shown in Fig- 
ure[5fb). About 85% of sgzKs in our sample {K^s, < 22) 
have SFR > 7OM yr~ 1 , and the median SFR is about 
184 Mq yr -1 , which is similar to the values in Daddi 
et al. (2004) (with a typical SFR of ~ 200 M e yr" 1 ) and 
in Kong et al. (2006) (with a median SFR of - 190 
Moyr- 1 ). 

4.3. The SFR—M* correlation 

The existence of a correlation between star formation 
and stellar mass of galaxies at different redshifts (0 < z < 
2) has been reported recently (Daddi et al. 2007; Elbaz 
et al. 2007; Pannella et al. 2009). To gain further insights 
into the nature of star formation at 1.4 < z < 3.0, we 
show the SFR as a function of M* in the four redshift 
bins in Figure [HI To examine in detail the stellar mass 
dependence of the SFR, we plot a good-fit line (dashed 
line) for galaxies with M* > 1O 1O M0 in each redshift 
panel. From this figure, we find a positive correlation 
between the SFR and stellar mass of sgzKs in the all 
four redshift bins. 

In Figure [Sfd), we overplot the correlation between the 
SFR and stellar mass for z ~ 0.1 galaxies in the Sloan 
Digital Sky Survey (SDSS; Brinchmann et al. 2004), z ~ 
1 and z ~ 2 galaxies in the GOODS (Elbaz et al. 2007; 
Daddi et al. 2007), and z ~ and z ~ 1 from Millennium 
model (Kitzbichler & White 2007). Figure EJd) shows 
the slopes between the SFR and stellar mass of galaxies 
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Fig. 4. — X-band differential number counts for star-forming galaxies (sgzKs) and passively evolving galaxies (pgzKs) at z ~ 2. For 
comparison, we have overplotted the number counts for z ~ 2 galaxies from the literature and the predictions of model of Kitzbichler & 
White (dot-dashed line). Left: sgzKs. Right: pgzKs. 



O 

o h 



O 
• O |- 

E ^ 
3 



o 
o h 



-i — i — i — r- 



a) 



T 1 1 T- 



T 1 1 1 - 



T 1 1 1 - 



jzm 



o 

O |- 



c O 
£= o 

3 CM 



O 

o 



b) 



T 



T 1 1 T- 



T 



t — i — i — r 



T 



T 1 1 1 - 



10 



11 



12 



13 



log(MyM ) 



log SFR (M yr" 1 ) 



Fig. 5. — Stellar mass and star formation histogram of sgzKs in the EGS field. The left panel (a) is the plot for stellar mass, and the 
right panel (b) is the plot for star formation rate. 

at different redshifts are very similar. On average, the 
SFR of galaxies on the SFR— M* diagram increases with 
redshift at a fixed M*, by a factor about 6 from z ~ 
to z ~ 1, and a factor about 20 from z ~ to z ~ 2. At 
fixed stellar mass, star-forming galaxies were much more 
active on average in the past. This is most likely due to 
a larger abundance of gas, depleted with passing of time 
(Daddi et al. 2007). 

The Millennium model overlaps very well with the 
SDSS at z ~ 0.1 and predicts a slope very similar to 
the observed trend at different redshift bins. For z <~ 1, 
however, we find that at fixed stellar masses the model 



predicts an increase of the SFR much below the observed 
trend. It seems that a major change required for mod- 
els would be to increase the star formation efficiency 
(and thus the typical SFR) at all masses for star-forming 
galaxies at redshifts z > 1. 



4.4. The sSFR-M* correlation 

sSFR is defined as the ratio of the current SFR to 
the current stellar mass, sSFR=SFR/Af*. Thus, higher 
values of the sSFR indicate that a larger fraction of stars 
were formed recently, and it can be used to characterize 
the star formation history of galaxies. 
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Fig. 6. — Relationship between stellar mass with SFR, for z ~ 2 star-forming galaxies in our sample for each redshift bin. Dashed 
lines and the equations ( here Y = log SFR, X = log (M* / Mq)) in each panel show the linear fit for galaxies with M* > 10 10 Mq. The 
dot-dashed lines are the z = 2 and 2 = 1 correlations from Daddi et al. (2007) and Elbaz et al. (2007). The solid lines show a prediction 
for 2 = and 2 = 1 from the Millennium simulations (Kitzbichler & White 2007). The dashed line is the z = 0.1 correlation from SDSS 
(Brinchmann et al. 2004). 

the lowest sSFR at redshift z < 2; they formed their 
stars earlier and more rapidly than their low-mass coun- 
terparts. This is in agreement with the downsizing sce- 
nario. 



Figure [7] shows the sSFR as a function of stellar mass 
for the sgzKs in our sample, with the same four red- 
shift bins as in Figure [5J The dashed lines in each panel 
represent a good fit between the sSFR and Af* for galax- 
ies with > 10 10 M Q . In Figure E[d), we overplot 
the correlation between the sSFR and stellar mass for 
0.2 < z < 0.5, 0.5 < z < 0.75, and 0.75 < z < 1.0 
galaxies in Brinchmann & Ellis (2000). 

From this figure, we find that the general shape of 
sSFR— M* correlation is in good agreement with a simi- 
lar study for z < 1.0 galaxies (Brinchmann & Ellis 2000), 
and the line of maximum sSFR runs parallel to lines of 
constant SFR. Furthermore, this sSFR is generally in- 
creasing with increasing redshift for all stellar masses, 
suggesting that star formation contributes more to the 
growth of low-mass galaxies than to the growth of high- 
mass galaxies. Finally, the most massive galaxies have 



5. MORPHOLOGY 

The morphology of a galaxy reflects its dynamical his- 
tory and evolution. For instance, different Hubble types 
are associated with different star formation histories and 
different patterns of motion of stars and gas. Further- 
more, galaxy morphology correlates with a range of phys- 
ical properties in galaxies, such as stellar mass, lumi- 
nosity, size, velocity dispersion, concentration, and, par- 
ticularly, color; this suggests that morphology is crucial 
in our understanding of the formation and evolution of 
galaxies (Fang et al. 2009; Kong et al. 2009). 

With the advent of the new WFC3 on board HST, 
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Fig. 7. — Specific star formation rate (sSFR) as a function of stellar mass for sgzKs in the EGS. The dashed lines correspond to the best 
fit to the four redshift ranges indicated. The three dashed lines from Brinchmann & Ellis (2000) were overplotted in the bottom-left panel 
(d). The sample of sgzKs with open circles from the overlapping area between NMBS and AEGIS. 



with its vastly improved sensitivity and resolution com- 
pared to NICMOS, it has become possible to analyze 
the rest-frame optical structure of z ~ 2 galaxies with 
an unprecedented level of detail. To describe clearly the 
morphologies of both sgzKs and pgzKs, we have per- 
formed nonparametric measures of galaxy morphology by 
using data from HST WFC3/H(F160W) imaging (Gro- 
gin et al. 2011; Koekemoer et al. 2011), such as Gini co- 
efficient (the relative distribution of the galaxy pixel flux 
values, or G) and M20 (the second-order moment of the 
brightest 20% of the galaxy's flux). Early- type galaxies 
have higher G and lower M 2 o, while late- type galaxies 
have lower G and higher M2o- Combining G and M20 
permits classifying different galaxy population effectively 
(Lotz et al. 2006; Kong et al. 2009). 

In Figure [8j the blue squares and red circles repre- 
sent the morphological distribution of 173 sgzKs and 50 
pgzKs (these sources have H-band counterparts) in the 



rest-frame wavelength ~ 5300A, respectively. The distri- 
bution of gzKs is very similar to that of local galaxies (E- 
Sd) in Lotz et al. 2004, with pgzKs showing high G and 
low M20, and sgzKs with lower G and higher M20 values. 
This figure clearly illustrates the morphological variety 
present at z ~ 2. The median (G, M20) values for sgzKs 
are (0.49, -1.48) in the rest-frame optical band, while 
pgzKs are (0.64, -1.74). A sample of 73 local ULIRGs 
(z < 0.2) is also shown in this diagram (Lotz et al. 2004), 
indicated as empty stars, using data from HST WFPC 
2/F814W imaging (Borne et al.2000). The solid dia- 
monds and triangles in Figure [8] represent the starburst- 
dominated sources at z ~ 2 from Bussmann et al. (2011) 
and Forster Schreiber et al. (2011), respectively; their G 
and M20 values were derived from NICMOS/F160W im- 
ages, corresponding to the rest-frame optical wavelength. 
The dashed line (G = 0.4A/20+0-9) in this figure is drawn 
qualitatively at rest-frame optical band (~ 5300A; Buss- 
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Fig. 9. — HST/WFC3 J- and //-band images of six sgzKs in the 
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M20 > —1-10. Values of G and M20 are shown in each //-band 
image. All images are in negative gray scale and are 4x4 arcsec 2 
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M 20 

Fig. 8. — M20 vs. Gini coefficient (G) for sgzKs, pgzKs, z ~ 2 
ULIRGs and local ULlRGs (empty stars; Lotz et al. 2004). The 
dashed line is defined as G = OAM20 + 0.9 (Bussmann et al. 2011). 
Filled blue squares and red circles represent the morphological 
distribution of 173 sgzKs and 50 pgzKs in our sample, respec- 
tively, which are matched in //ST/WFC3 //-band. The starburst- 
dominated sources at z ~ 2 from Forster Schreiber et al. (2011) and 
Bussmann et al. (2011) are shown as filled triangles and diamonds, 
respectively, with G and M 20 derived from HST NICMOS/F160W 
images. 

mann et al. 2011). Galaxies with G < 0.4M 20 + 0.9 have 
diffuse structures or multiple bright nuclei in appearance. 
Objects with G > 0.4Af 2 o + 0.9 and -1.6 < M 20 can be 
irregular, often with a bright nucleus with tidal tail or 
fainter knots. As for sources with G > OAM20 + 0.9 
and M20 < —1.6, they are relatively smooth with a sin- 
gle nucleus. From Figure [5J we find that star-forming 
galaxies at high redshift not only have diffuse structures, 
similar to starburst-dominated ULIRGs at z ~ 2 with 
lower G and higher M20 values, but also have single- 
object morphologies (higher G and lower M20). More- 
over, the distribution of M20 of sgzKs is mainly located in 
the range from -1.0 to -1.8; these features of morpholo- 
gies are consistent with expectations from simulations 
of major mergers during the beginning and end stages, 
respectively, of the phase of the merger when the SFR 
peaks and begins to turn over (Lotz et al. 2008). 

Figures 1 and OH show the HST/WFC3 J- and H- 
band stamp images for a subset of six sgzKs and six 
pgzKs, respectively. Their rest-frame optical morpholo- 
gies for sgzKs are very diversified including string-like, 
extended/diffused, and even early-type spiral morpholo- 
gies, while pgzKs are relatively smooth and compact, im- 
plying that there are different formation process for these 
galaxies. Therefore, we conclude that the population of 
Hubble sequence galaxies roughly matches that of the 
peculiars sometime between z = 1.5-2. 

6. FRACTION OF AGNS 

Mid-IR photometry has been proven to be a robust 
and efficient tool to select AGNs without prior infor- 
mation (e.g., Lacy et al. 2004; Stern et al. 2005; Donley 
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et al. 2008), as their properties at these wavelengths are 
typically very different from those of stars and galaxies. 
There have been several studies of AGN selection using 
mid-IR color or other IR properties. In this section, two 
different AGN selection methods are employed, to iden- 
tify AGNs in our gzK sample. The first one is based 
on the IRAC band color criteria of Stern et al. (2005). 
The second method is based on Mid-IR spectral index 
(Barmby et al. 2006; Park et al. 2010). 

6.1. Mid-IR colors 

Figure QT] shows the mid-IR color space of the EGS 
survey along with the Stern et al. (2005) AGN selection 
criterion (green dot-dashed lines). Small dots represent 
the distribution of Mid-IR color for sources detected in 
all four IRAC channels in the AEGIS field; they have 
very similar [3.6] - [4.5] colors, as this is dominated by 
their stellar emission. Six hundred and twelve sgzKs in 
our sample are detected in all the four IRAC channels 
and are plotted as blue open circles. One hundred and 
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Fig. 11. — IRAC band color-color diagram of galaxies in the 
EGS field. Dot-dashed lines indicate the boundaries of the AGN 
selection wedges from Stern et al. (2005). Small points represent 
all galaxies, sgzKs are denoted by blue open circles, and pgzKs are 
denoted by red open circles. 

twenty-two of them can be selected as AGN candidates 
by the mid-IR selection criterion of Stern et al. (2005). 
Two hundred and sixteen pgzKs in our sample are de- 
tected in all the four IRAC channels and plotted as red 
open circles. Thirty-five of them are AGN candidates. 
Therefore, using the mid-IR AGN selection criterion, 
~ 20% of sgzKs and ~ 16% of pgzKs are classified as 
AGNs, respectively. 

The AGN selection criterion from Stern et al. (2005) 
were based on the spectroscopic sample of the AGN and 
Galaxy Evolution Survey, most of which have redshift 
z < 0.6. Their mid-IR color criterion is reliable for 
classifying AGNs and galaxies at low redshift. In con- 
trast, observations and templates suggest that a high de- 
gree of stellar contamination is unavoidable when the 
color selection technique is applied to deeper samples 
(Barmby et al. 2006; Donley et al. 2008; Georgantopou- 
los et al. 2008; Assef et al. 2010). The redshift of sgzKs 
and pgzKs in our sample is z ~ 2; with much deeper 
IRAC data, we will classify AGNs from normal galax- 
ies with the mid-IR spectral index method in the next 
subsection. 

6.2. Mid-IR spectral index 

The slope of a galaxy SED in the infrared can be 
characterized by a power-law behavior of flux density 
with frequency /„ oc v a . AGN SEDs often follow a 
negative-sloping (red) power-law, which may arise from 
either thermal or non-thermal emission originating near 
the central region. In contrast, stellar-dominated sources 
at redshift z ^ 2 generally exhibit positive (blue) IRAC 
power law emission. For this reason, mid-IR power-law a 
of galaxies have been proposed and used as classification 
criteria for AGNs and normal galaxies, with a small level 
of galaxy contamination. 

In this paper, the IRAC fluxes for each gzKs, covering 
the 3.6-8.0 range, are fitted with a power law, a. 



The a value is only accepted if the \ 2 probability fit is 
P x 2 > 0.1. Galaxies fitted well with a power-law SED are 
further classified as red or blue. A limit of a iS — 0.5, 
following Donley et al. (2008) and Park et al. (2010), is 
chosen to classify gzKs as AGNs. The criterion reveals 
that the fraction of AGNs in our sample is ~ 10%. In 
the meantime, we also check the X-ray luminosities of 
these AGN candidates (if they have X-ray counterparts) , 
and the result shows that most of these objects have 
Ln. 5-1 n koV > 10 41 ' 



erg s " . 

7. SUMMARY 



In this paper, we have described the construction of 
a sample of gzKs (correspond to the BzKs of Daddi 
et al. 2004) within the AEGIS. With this field's unique 
combination of area and depth at almost every waveband 
observable, we analyze the photometric redshifts, num- 
ber counts, SFRs, sSFRs, stellar masses, morphologies, 
and the fraction of AGNs of gzKs in the AEGIS. Our 
main conclusions are as follows. 

1. We identify a sample of 2031 gzKs including 1609 
star-forming galaxies, with gzK = (z — A')ab — 1-4(5 — 
z)ab > 0.2 and 422 passively evolving galaxies with 
gzK < 0.2 and (z - K) AB > 2.7, to a limit K AB < 22. 
The surface density of sgzKs+pgzKs is ~ 1.28 arcmin -2 . 
Of the sgzK galaxies, 83.2% lie in the range 1.4 < z < 
2.5, while 8.5% are at z < 1.4 and 8.3% at z > 2.5. For 
pgzK galaxies, there is a considerable fraction (12.8%) of 
objects at low redshift (z < 1.0). The z p median values 
of both sgzKs and pgzKs are 1.8 and 1.5, respectively. 

2. For a sample selected at Kab < 22.0, 2989 spec- 
troscopic redshift of field galaxies (including 35 LBGs at 
z ~ 3 and 14 ULIRGs at z ~ 2), we reach an average red- 
shift accuracy of 6z/(l + z s ) = —0.014 with the normal- 
ized median absolute deviation (ctnmad) of 0.032. Only 
4.3% galaxies have (z p — z s ) / (1 + z s ) > 0.1. We also com- 
pare our photometric redshifts with NMBS photometric 
redshifts in the AEGIS field; they are in good agreement. 

3. Similar to the findings in other literatures, our 
counts of passively evolving galaxies turn over at K^a ~ 
21.0 and both the number of faint and bright sources 
(including sgzKs and pgzKs) in our catalogs exceed the 
predictions of a recent semi-analytic model of galaxy for- 
mation, implying that a more successful model must have 
to explain this diversity in the future. 

4. Based on the reddening-corrected UV luminosi- 
ties and SED fitting of the multi-band photometry, we 
find that gzA-selected star-forming galaxies have a me- 
dian SFR and stellar mass of ~ 184M Q yr _1 and ~ 
8.8 x 10 10 Mq, respectively. Moreover, we also find that 
the SFR and sSFR of sgzKs increases with redshift at all 
masses, suggesting that star- forming galaxies were much 
more active on average in the past. Comparing to the 
predictions of the Millennium model, there is a major 
change required for models in order to increase the star 
formation efficiency at all masses for star-forming galax- 
ies at redshifts z > 1. 

5. We also study the morphological properties of gzKs 
in our sample, employing data from HST WFC3/F160W 
imaging within the EGS field. We derived the median (G, 
M20) values for sgzKs and pgzKs are (0.49, -1.48) and 
(0.64, -1.74) in the rest-frame optical band, respectively. 
Moreover, we find that morphologies of z ~ 2 galaxies 
in our sample are complex and varied: from compact, 
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apparently early-type galaxies to large star-forming sys- 
tems superficially similar to pure disk or irregular galax- 
ies, implying that there are different formation process 
for these galaxies. We conclude that the population of 
Hubble sequence galaxies roughly matches that of the 
peculiars sometime between z =1.5-2. 

6. For the fraction of AGNs in the gzKs, we compared 
the power-law selection technique with mid-infrared color 
selection, and we find that 82 sources with a < —0.5, 
corresponding to the fraction of AGNs, is ~ 10% in our 
sample of 828 gzK galaxies with four IRAC bands. In 
the meantime, we also check the X-ray luminosities of 
these sources (if they have X-ray counterparts); the re- 
sult shows that most of these objects have L 5-10 kcV > 
_i 

erg s . 
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